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Abstract 

We will consider two special families of polynomial perturbations of the linear center. For the resulting perturbed 
systems, which are generalized Lienard systems, we provide the exact upper bound for the number of limit cycles 
that bifurcate from the periodic orbits of the linear center. 



Keywords: 

limit cycle, Lienard system, periodic orbit 



1. Introduction and statement of the results 

The bifurcation of limit cycles by perturbing a planar system which has a continuous family of cycles, i.e. periodic 
orbits, has been an intensively studied phenomenon; see for instance [3] and references therein. The simplest planar 
system having a continuous family of cycles is the linear center, and a special family of its perturbations is given by 
the generalized polynomial Lienard systems: 

x = y + Y, e 'W' y = Z eigi(x) ' 

i=l ;=o 

where jj. e N, v e N U {0}, go(x) = -x, gi(x) and Fj(x) are polynomials for i > 1, and s is a small parameter. 

The classical and generalized Lienard systems appear very often in several branches of science and engineering, 
as biology, chemistry, mechanics, electronics, etc. see for instance lil2ll and references therein. In particular Lienard 
systems are frequent specially in physiological processes, see for instance [5]. In addition, the family of generalized 
polynomial Lienard systems is one of the most considered families in the study of limit cycles, see 1111. 

We assume that F^x) £ 0, g v (x) £ 0, m = maxi<,< / ,{deg,F,'(x)}, and n = maxi<,< v {deg,g,(x)}. For a small enough 
s, let "Ky(m, ri) be the maximum number of limit cycles of ( fQ ) that bifurcate from cycles of the linear center (lo), i.e. 
the maximum number of medium amplitude limit cycles which can bifurcate from (lo) under the perturbation ( fQ l. If 
v = 0, then Ji^(m, n) does not depend on n; hence we only write The main problem concerning "Ky{m, n) is 

finding its exact value. 

We know from [9] that "K^im) > [(m - l)/2], where [■] denotes the integer part function. Moreover, by following 
J3, Theorem 3.1] we can prove that ^(m) = [(m - l)/2] for /j. > 1; Theorem Q] (below) is a generalization of this 
result. Also, we know from that JC\(m,ri) > [(m - l)/2], "K\{m,n) > max{[(m - l)/2] , [m/2] + [n/2] - 1}, and 
J£o(w, n) > [(m + ri)/2] - 1. However, the exact values of "K \{m, n), r Ki(m, «), and l Kl(m, n) were not reported there. 
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In this paper we give the exact value of Jfy(m, n) for two subfamilies of ( fQ ). More precisely, we will give the 
exact value of !Ky(m,«) and 5iy{m,n), where 3{.y(m,n) is the value of r Ky{m,ri) by assuming that gj(x) is odd for 
1 < i < v, and 3-Cy(m,n) is the value of r Ky(m,ri) by assuming that Fi(x) is even for fio < i < p, where fio with 
1 < IM) < fiis the smallest integer such that F^x) £ 0. Of course, if fio = fi, then "Kyim, n) = 3{.y(m, n). 

Our main result is the following: 

Theorem 1. (a) "K^m, n) = fel. (Z?) SC^m, n) w either ifm is odd or [f ] + [f ] - 1 i/m w even. 

The assumptions on and Fj(x) in definitions of l Ky{m, n) and J£y(m, n), respectively, are necessary. Other- 
wise, we can construct systems ( fQ l having more medium amplitude limit cycles, see Remark 1 in Section[3] 

Theorem[Tj£>) is a generalization of Theorem 1 . 1 in II 1311. w here the case yi-v-X was considered. We note that 
in such a case 5i\(m,n) = Jf[(m,n). Hence Theorem 1(b) ([13, Theorem 1.1]) gives the exact value of "K\(m,ri). 

The proof of Theorem[T]is based on computing the maximum number of isolated zeros of the first non-vanishing 
Poincare-Pontryagin-Melnikov function of the displacement function of fQ i, by taking into account the restrictions: 
gi(x) odd for 1 < i < v and Fj(x) even for (/$ < i < fi, respectively. 

The paper is organized as follows. In Section [2] we recall the definition of the displacement function of (fO, as 
well as the algorithm to compute the Poincare-Pontryagin-Melnikov functions. Preliminary results that allow us to 
provide elementary proofs of the main results are given in Section[3] Finally, in Section|4]we will prove TheoremQ] 

2. Poincare-Pontryagin-Melnikov functions 

The linear center (lo) is the Hamiltonian system associated to the polynomial H - (x 2 + y 2 )/2; hence its cycles 
are the circles y c — {H - c — 0} with c > 0. By using c as a parameter, the first return map of ( fQ l can be expressed 
in terms of s and c: [P(e, c). Therefore the corresponding displacement function Lie, c) = T(e, c) - c is analytic for 
small enough s and can be written as the power series in s 

L(e, c) = eLi(c) + e 2 L 2 (c) + 0(s 3 ), (2) 

where L,(c) with i > 1 is the Poincare-Pontryagin-Melnikov function of order i, which is defined for c > 0. 

Let Lic(c) with k > 1 be the first non-vanishing coefficient in (0. The zeros of L^ic) are important in the study 
of medium amplitude limit cycles of fQ i because of the Poincare-Pontryagin-Andronov criterion: The maximum 
number of isolated zeros, counting multiplicities, of Lt(c) is an upper bound for 'Ky(m, n). Furthermore each simple 
zero co of L^ic) corresponds to one and only one limit cycle of fQ ) with s small enough bifurcating from the cycle 

We know from [6] that Lk{c) has at most [£(max{n,m) - l)/2] positive zeros, counting multiplicities. However, 
this result does not give the value of "Kyim, n) because the upper bound for k depending on fi, v, m, and n is unknown. 
Now, we will recall the algorithm to compute the functions L,(c). System ( fQ ) can be written as 

x=y, y = -x + e{g\ (x) + f (x)y) + e 2 (g 2 (x) + fi(x)y) + ■■■ , 

where fix) = F'Ax), or equivalently as 

dH - sa>i - e 2 L02 - • • ■ = with oj, = (gi(x) + f(x)y) dx and u>j = for i > maxf/v, v}. (3 £ ) 

As we know, L\ (c) is given by the classical Poincare-Pontryagin formula L\ (c) = f io\. The result for computing 
the higher order Poincare-Pontryagin-Melnikov functions is the following: 
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Theorem 2. (Yakovenko-Francoise-Iliev algorithm H], 10]). > 2 and L x {c) = ••• = L*_i(c) = 0, ffte« 
f/iere are polynomials q\,..., qt-i and Qi,. . ., Qk-i such that Q[ = dQ\ + q\dH, . . . , Qt_i = dQk-i + qt-idH, and 



L k (c)= f Q*, 

£2i = a>i, awaf D.[ — (±>i + ^ g/W/ w/f/i /, y > 1 /or 2 < I < k. 

<+H 

The proof of this result easily follows from the Poincare-Pontryagin formula, and the Ilyashenko-Gavrilov theo- 
rem (Q, [Q]]): If JT, <w = for all c > 0, then <w = dQ + gafi/, where 2 and q are polynomials, and by applying an 
induction argument. For a detailed proof, see for instance @], |0]. 

To simplify the computation of the Poincare-Pontryagin-Melnikov functions, we will give some properties of w,-. 



3. Preliminary results 

For computing L^(c) for fQ we will use the following two elementary lemmas whose proof is omitted. 
Lemma 3. Let P be a polynomial in the ring R |x 2 , //]. We define deg 2 P to be the degree of P in R[x 2 , H]. 

(a) For i, j > there are homogeneous polynomials Q,j, q,j € K |x 2 , i/J vv/f/i deg 2 <2y = i+j anddeg 2 q,j — i+ j—l, 
such that H l x 2j dx = d (xgiy) + (xqijj dH or H'x 2j+l dx = d (x 2 g<y) + (x 2 <?,y) dH. Ifi = 0, then qy = 0. 

(b) For i, j > there are homogeneous polynomials Qy, qij € R[x 2 , H] with deg 2 Qij — i + j+l and deg 2 qij — i + j, 
such that H l x 2 i +i ydx — d (yQijj + (yqij) dH. 



(c) For i, j > we have 



f //'x ; 



.'+./ 



2J(2j + l)\ j+l j 

Lemma 4. Ifaj e A:- {(xA + xyB)dx\ A, Z? e R^x 2 ,//]} anofg eS:= {x 2 gi +y^2 1 q\,qi e R [* 2 ,//]}, then qu> e A 

The next two results are straightforward consequences of these two previous lemmas. 
Corollary 5. Ifaj e A, then J a> = 0, d — dQ + qdH with q e §, and qa> e A. 



t d 



z 

Vr=0 



2(r+l) 



r+1 /2 r (2r+l) 



Corollary 6. IfP{x 2 ) = V /? f x 2 '' e R [x 2 ], then J p{x 2 )ydx = -nc 

The following two lemmas will be important in the proof of TheoremQ] 
Lemma 7. Suppose Theorem\2\ IfCli e A for 1 < I < k — \ , then a>i e A for 1 < I < k — 1, and Lk(c) — J" w^. 

Proof. We proceed by induction on £. If A; = 2, then fii = Wi e A Hence = af(2i + q\dH, q\cj\ 6 A, and 
JT q'iWi = by Corollary|5] Since L2(c) — Q2 = J u>2 + f y qi^i because of Theorem|2] L,2(c) = JT w 2 . 

We assume that the lemma is true for k—1, and we will prove it for k. By assumption, Sli e A for 1 < I < k — 1. 
Then, by Corollary [5] Q; = + <37<£ff with <?/ e S for all 1 < / < it — 1. In addition, by the induction hypothesis, 
aii e A for 1 < I < k - 2. Thus, Qt_; := 2i+j=*-i ^i<w/ with i, j > 1 is an element of ./I following Lemma[U Since 
<x>k-\ = Q/t-i - Q-k-u (jJk-\ 6 A Hence it is clear that := Yit+j=k Qi^j w i m ! > 7 ^ 1 i s an element of A, which implies 
that Q.k = by Corollary |5] Finally, from Theorem[2] we have L^c) = j£ Wyt + J^, Therefore La(c) = u>k. □ 



Before announce next lemma, we note that each polynomial h(x) = zZ'"=o a r x r of degree m — 1 can be written as 



[¥] _ [¥] 

h(x) = /z(x 2 ) + xh(x 2 ^, where h(x 2 } = ^ az r +\X 2r , and h(x 2 ^ = ^ ci2r+2X 2r . 

r=0 r=0 

Lemma 8. Lef w = (g(x) + f(x)y) dx, where f(x) = Yj7=o a r xr an ^ 8( x ) - Z"=o b s x s . 

\ r ) 
(b) If J co = 0, then co — dQ + (yq)dH with ^el |x 2 , i/J of degree deg 2 q = [(m - 2)/2], anc/ 



(4) 



(^2. ; )(a2r+2) 

2- !+r (2s + l/ 



( c ) (y?) w = if and only if q = 0, or g (x 2 ) = 0. 

Proof, (a). By (a) and (£>) of Lemma[3] a) = f /(x 2 )yofx. Finally, the statement follows from Corollary|6] 

(b). If jT co = 0, then /(x 2 ) s by statement (a). This property implies that w = g(x)dx + xf(x 2 ^ydz 
d{J g(x)dx^j + z[=o ^ fl2r+2X 2r+1 y^ x by ©■ From LemmaFJJZ?) we obtain x 2r+x ydx = c/(y(2r) + (yq r )dH, thus 



oj = d 



I 



g(x)ofx + y 



2^ a 2r+lQ r + y °2r+2q r 



dH = dQ + (yq) dH, 



where Q r , q r , q e Rjx 2 , //J are homogeneous and deg 2 q — F 2 ^]- Moreover, a simple computation shows that 



9r 



r+l\(r+l—i\, ; / , \i 
— \(2H) r -' (x 2 -2H) 



i I \ 2i 



(5) 



As (yq) co = qg(x 2 ^jydx + qgix^jxydx + qf(x 2 ^jxy 2 dx and qf[x L ^xy 2 dx = (7/(x 2 )x (2// - x 2 ) c/x, it follows that 
(yq) to = qg{x 2 )jydx + dQi + qidH because of statements (a) and (b) of Lemma|3] Hence we obtain 



J, r. Vm-2)/2] L „,z.j 

(yg) co = I ^(x 2 )yt/x = I ^ fl2r+2? r ^ fo 2. s x 2s 



yc/x. 



By using expression (0 of q r , a straightforward computation, and Lemma [3jc) we obtain the formula given in the 
statement. Finally, statement (c) follows from the formula given in statement (b). □ 

Remark 1. System ( fQ > with fi — v = 1, Fi(x) = -x 2 , and g\(x) — 1 - x 2 does not satisfy the hypothesis in definition 
of 5iy(m, n) because g\(x) is not an odd function. Here m — n — 2 and from TheoremQJa) it follows that !K}(2, 2) = 0; 
however, for e small enough, this system has one medium amplitude limit cycle. Indeed, we need only to prove that 



the first non-vanishing coefficient of the displacement function (0, associated to the system, has a simple positive 
zero. The system can be written in the form ( |3~e| as dH - e<±> — with a> = (1 - x 2 - 2xy)dx. By Lemma Oa), 
L\(c) = 0, and by Theorem[2]and LemmaOb), L2(c) = -nc(4 - 2c). Now, system ( fQ ) with fi — v — 2, F\{x) — -3x 2 , 
Fi(x) - -2.x 3 , g\(x) — x 2 + x 3 , and gz{x) = (-5 + 25x 2 ) /6 does not satisfy the hypothesis in definition of ?Cy(m, n) 
because /^(x) is not an even function. In this case m = n = 3 and by Theorem[TJ^), 5^(3, 3) = 1; however, for s 
small enough, the resulting system has two medium amplitude limit cycles. Indeed, following previous ideas, and 
using Theorem[2]and Lemma|8]it is easy to see that L\{c) = 0, L^ic) = 0, and Lj,(c) = -nc(c — l)(c - 2). 

4. Proof of the main results 

We can assume, after a linear change of variables if necessary, that F,(0) = for all 1 < i < fi. Suppose that 
Fi(x) = Z^Li( a i(r-r>/r)x r and g t (x) = E"=o*w^- Thus ' M x ) = F 'i( x ) = IT=o a irX r and g t (x) can be written as 
fi(x) = fi{^) + x fi( x2 ) an d gi( x ) = g ;(* 2 ) + xg i(* 2 ), respectively, according to ©. 

Proof of Theorem\T\ (a). By hypothesis, is odd for 1 < i < v, which means that gi(x) = xg ,(x 2 ) for 1 < i < v. 
Let Lic(c) be the first non-vanishing Poincare-Pontryagin-Melnikov function in (0. If k = 1, then the theorem is 
true. Indeed, we have L\(c) — coi — J xg {(jrjdx + f ^(x^ydx + J f\[x 2 ^xydx, and as xgi\x 2 jdx = 0, and 

J y f\{x 2 ^xydx = by Corollary [5] we obtain L\(c) = f\[x 2 ^ydx. From (O we have deg 2 /i(x 2 ) = [(m - l)/2]; 
hence L\(c) has at most [(to - l)/2] positive zeros because of Corollary [6] Moreover, we can choose suitable coef- 
ficients of F\{x) in such a way that L\(c) has exactly [(to - l)/2] simple positive zeros. Therefore, by applying the 
Poincare-Pontryagin-Andronov criterion it follows that !Ky(m, n) - [(to - l)/2]. 

Suppose then that k > 2 and we are therefore in the hypothesis of Theorem [2] If Q/ £ A for 1 < / < k — 1, 
then Lk(c) = JT o>* by Lemma [7] and by applying the same idea as in previous paragraph, we obtain $Cy(m,n) = 
[(to - l)/2]. Accordingly, it remains to prove that Q/ E A for 1 < I < k — 1. 

We proceed by induction on k. If k — 2, then Li(c) = 0, which implies that Qi = (xgi(x 2 ) + xyf\\or\\dx E A. 
We now assume that the assertion is true for k - 2, and we will prove it for k - 1. By induction hypothesis, Q, e j4 
for 1 < i < k - 2, which implies that Q, = dQi + q/dH with o,- e S for 1 < i < k — 2 by Corollary [5] Furthermore, 
by Lemma|7l o> ; - E ./I for 1 < j < k - 2. Hence Q^-i : = Y^i+j=k-\ Qi^j w i m 1 < i, J < A — 2 is an element of 
^4 because of Lemma |U Since f2^_i = cj^-i + Q^_i, Lj:-i(c) = J Q^_i = J u>k-i = fk-\(x 2 ^ydx = 0, whence 

(*>k-\ = (xgk-iyx 2 ) + ■Q'A-i^ 2 )) dx E A Therefore Oj-i e A, which completes the proof of statement (a). 

(b). First, we will note two properties concerning a>j and j£ u)i which we will use along the proof. Second, we 
will split the proof into two cases: m odd and m even. 

For 1 < i < hq the 1-form w, = gj(x)dx is exact, that is, w, = dQi + qjdH with q, = 0. Hence, by Theorem [2] 
Q, = a>i and L,(c) = f Q; = for 1 < / < //q, and L^c) = = w^. On the other hand, since F t {x) 

is even for fio < i < /i, f(x) = xfi^x 2 ^ for fio < i < yL. Thus, w,- = d{J gj(x)dx^ + xf(x 2 ^jydx for i > fiQ, and as 
x 2r+l ydx — d {yQo r ^ + (yqo r ) dH because of Lemma[3f^), we conclude that a>i = d \Qj} + (yq,) dH; of course qi = 
for i > yu. Therefore J a); = for all i > //o- 

Case m odd. If ;« is odd, then deg/^ (x) = m because Fj(x) is an even polynomial for fto < i < fi. Since 
■^/ioW = ffo^x) - fn{ x2 ) + Jt: ^o(- ;ic2 ) nas an even degree, ^ (x 2 ) ^ 0. Hence, from Lemma [HJa) it follows that 
Lfi (c) = JT cj^ = J /^(x^ydx ^ 0, and it has at most [(m - l)/2] positive zeros, counting multiplicities; moreover, 
we can choose suitable coefficients of F m (x) in such a way that L^ic) has exactly [(to - l)/2] simple positive zeros. 
Therefore by the Poincare-Pontryagin-Andronov criterion, l Ky{m,ri) = [(to - l)/2]. 
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Case m even. Let Lk(c) be the first non-vanishing Poincare-Pontryagin-Melnikov function of (0. If k — yuo, 
then Lfi (c) has at most [(m - l)/2] positive zeros, counting multiplicities, because of Lemmata). Since m is even, 
[(m - l)/2] < [m/2] + [n/2] - 1. Hence L Mo (c) has at most [m/2] + [n/2] - 1 positive zeros, counting multiplicities. 

We claim that if k > fio + 1, then tot,..., Uk-i-^ £ A, £2, = dQi + qidH with qi e S for juo < i < A — 1, and 
Lk(c) = J (y<lti ) <^k-iM, = J^, ^g*-/^* 2 ))'^*- By assuming that this assertion is true and by applying Lemma ^b) 
we conclude that Lk(c) has at most [m/2] + [n/2] - 1 positive zeros, counting multiplicities; moreover, we can choose 
suitable coefficients of q^ and g k-^yx 2 ) in such a way that Lk(c) has exactly [m/2] + [n/2] - 1 simple positive zeros. 
Thus, by the Poincare-Pontryagin-Andronov criterion, 3{,y(m,n) = [m/2] + [n/2] - 1. Therefore, to finish the proof 
of statement (b) we need only to confirm the assertion, which we prove next by proceeding by induction on k. 

If k + 1 , then we will prove that £2^, = dQ n + q^dH with q m e S, and that L /Jo+ i(c) = q^g \{x 2 ^ydx. We 
know that £2 W = io^ , and from Lemma[8]^) it follows that £2^ = co^ = dQ^ + q^dH, where q^ = yq Mo SOeS. On 
the other hand, by Theorem[2] L w+ i(c) = jT£2 w+ i, where £2 w+ i =w / , 0+I +? I w (i0 + -'- + ^ _ I w 2 +? fl0 wi. Since =0 

for 1 < i < yu , £2^+1 = o> w+ i + Moreover, since f w wt] = 0, L^+ifc) = f = f <? w g i(jc 2 )yc/x. 

If = // + 2, then L^ 0+ i(c) = (y^„) wi = J 7 q^g^x^ydx = 0. Since ^ * 0, g i(jc 2 ) = by Lemma^c). 
This implies that £2i = cj\ e A, and by Corollary [5] £2; = t/Qi + q\dH with e S. Moreover, we know that 
W/j = dQ^ + q^dH with q m = yq m e §, and o> w+ i = + {yq^+^dH. Thus, £2 w+ i = + ^ Wi = 

+ Itia+idH with ^ 0+ i € S because of Corollary|5] On the other hand, from Theorem[2]we have 



«7y c . ^yc- ^y £ *Jyc 



As wi e .A and q^+i e S, then we have q^+iui e A. following Lemma|4]and JT g^+iwi = by Corollary|5] In 
addition, we know that q { = for 1 < i < [i and f w w+2 = 0. Hence L^ 0+2 (c) = (y^ ) w 2 = q^g 2 {x 2 pdx. 
We now assume that the assertion holds for k— 1, and we will prove it for £. By Theorem[2] Lk(c) = JT Q^, where 

Q.k= Uk+ q\U)k-\ + ■ ■ ■ + q^-iUk+i-^ + q^OJk-^ + q^+itok-i-no + • • • + qk-2^2 + 9k-i b) l- 

Since q t = for 1 < i < /u , £l k = a> k + q^k-^ + q^+i^k-i-^ +■■■ + qk-i^i + qk-i^i- 

On the other hand, from the induction hypothesis it follows that a>\,..., a>k-2-ti 6 A, Q, = dQ, + q/dH with qi e § 
for fio < i < k - 2, and L^_i(c) = J q M0 gk-\- m (xrjydx. Since Lk-i(c) = 0, gk-i-^x 1 ^ = because of LemmaUc), 

which implies that dJk-i-^ e A. Therefore, q^cjk-i-^ H H qk-^i + <lk-i^\ e ^1 by LemmaH] Moreover, we have 

(x>k-\ - d[Qk-\) + (yqk-\)dH, and by applying Corollary [5] we obtain 

CVi = cok-i + q^oJk-i-n + ■■■ + qk-^i + qk-iu\ = dQ k -\ + qk-\dH, with q k -\ e S. 

Hence g^+iwt-i-^o + ■ ■ ■ + qk-2^2 + qk-\U\ 6 -A by Lemma|4l In addition, cj* = t/^Qj;) + (yqk)dH. Thus, we obtain 

Part of the results of this work come from the author's postdoctoral stay at the Departament de Matematiques of 
the Universitat Autonoma de Barcelona. The author would like to thank the Centre de Recerca Matematica for their 
support and hospitality during the period in which this paper was written. 
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